Ionic criticality : an exactly soluble model 
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Gas- liquid criticality in ionic fluids is studied in exactly soluble spherical models that use interlaced 
sublattices to represent hard-core maWicomponent systems. Short range attractions in the uncharged 
fluid drive criticality but charged ions do not alter the universality class. Debye screening remains 
exponential at criticality in charge-symmetric 1:1 models. However, asymmetry couples charge and 
density fluctuations in a direct manner: the charge correlation length then diverges precisely as the 
density correlation length and the Stillinger-Lovett rule is violated at criticality. 

PACS numbers: 64.60.Fr, 61.20.Qg, 05.50.+q, 64.70.Fx 



The nature of gas-liquid (or, more generally fluid-fluid) 
criticality in systems in which long-range ionic interac- 
tions play a significant role has been a focus of attention 
since still-puzzling experiments questioned the appropri- 
ate universality class [ll. Beyond further experiments hi, 
numerous theoretical [2, Li U, U| and computational 
studies have been reported; however, basic questions re- 
main open. Certainly, the character of criticality depends 
on the range of the interactions: One expects an Ising 
critical point in a fluid with short-range couplings but 
mean-field behavior for interactions of sufficiently long- 
range. So, might the introduction of ions interacting via 
long-range Coulomb forces destroy an Ising-type critical 
point? Coulomb interactions are exponentially screened 
in a conducting classical fluid, as proved rigorously at low 
densities 0. Charge fluctuations in a fluid of S species 
of charges q a and valences z a = q a /q (where q is an ele- 
mentary charge), thus decay over a few Debye screening 
lengths fu, where 



* 'cT— 1 



(1) 



in which p = Pa is the overall density while the mole 
fractions are x a = p a / 'p. But, does exponential screening 
on this scale hold near and at criticality? 

Indeed, a major open issue is the behavior of the pair- 
wize charge correlations near the critical point, where the 
density fluctuations diverge strongly. With 



G aT {v) = (p (T (0)p T (r)> - p a pr 



(2) 



the correlation functions, G NN , G zz , and G NZ for the 
density, charge, and charge-density, are 0| 

G x ,(r;7» = ]T qf x qf* G„ T (r; T, p) , (3) 

z — '<T,T 

where X and Y may be N or Z while fl N = and i? z = 
1. Except at the critical point (T c ,p c ), itself, we may 
suppose that the corresponding structure factor, S NN , 
has the small k = |k| expansion 

S NN (V/S NN (0) = 1 + zZ p>1 (~) P C P P (T, P ) fc 2p . (4) 



Near criticality, S NN (0;T, p c ) diverges as 1/t 1 when t = 
(T — T c )/T c — > 0+, while the length £jv,oo characteriz- 
ing the exponential decay of G NN (r; T, p c ) diverges as 
i° N /t v (where short-range forces have been assumed). 
At criticality, density fluctuations are long-ranged and 
S NN (k;T c ,p c )~l/k 2 -*. 

By contrast, the charge structure factor should obey 

s zz (k) = o + e z .^ 2 Yl >9 H p ^y r ' p) k2p > ( 5 ) 

where the first vanishing term results from electroneu- 
trality reflecting the internal screening in an ionic fluid, 
while the Stillinger-Lovett sum-rule [a 



€z,i(T, p) — £ D (T, p) 



(6) 



characterizes the screening of external charges. Does this 
hold near and at criticality? Finally, we focus also on the 
charge correlation length £ z ,oo(T,p), that specifies the 
exponential decay of G zz (r) when r — > oo. How does 
iz,oa vary when £ N ,oo diverges near criticality? 

To obtain insight into these questions, we study mul- 
ticomponent lattice gas generalizations of the spherical 
model 0, EH E3| specifically designed to represent hard- 
core interactions and thus avoid the mutual "annihila- 
tion" of oppositely charged ions on the same site. This 
crucial feature, which (in contrast to ^jj) allows gas- 
liquid criticality to survive in the presence of Coulomb 
interactions, is accomplished by using a set of equivalent 
interlacing sublattices (with sites i at Rf with spacing a) , 
one for each of the a — 1,2, ... ,S distinct particle species 
[lflj . Thereby unlike charges cannot approach closer than 
an effective hard-core diameter a^: see, e.g., Fig. ^ 

To specify the models more fully, consider a multicom- 
ponent grand canonical lattice gas with site occupancy 
variables nf = 0, 1; this is equivalent to an Ising magnet 
with spins sf = 2nf — 1 = ±1, subject to fields h a (lin- 
early related to the chemical potentials, p a ). For attrac- 
tive interparticle potentials, the corresponding spin-spin 
couplings, Jcrr(R), at lattice separations R, are positive. 
We decompose these couplings as 



J ffT (R) = J° T (R) - ±<? CT ^ C (R) 



(7) 
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FIG. 1: Interlaced + and — sc sublattices with ao = |v3a. 



where J° T (H) represents short-range or, more generally, 
integrable interparticle interactions which we suppose 
suffice to drive gas-liquid criticality (even if the charges go- 
vanish). We take the <i-dimensional Coulomb potential, 
</? c (R) ~ 1/R d ~ 2 (d > 2), as the solution of an appropri- 
ate discrete Laplace equation 0] (with, for convenience, 
a uniform neutralizing background so that electroneutral- 
Ity, J2 a x <y ( l<y — 0) must be imposed). Fourier transforms, 
J aT Q&), are defined by summing over one sub-lattice, with 
Brillouin zone B, and it is useful to introduce 



AJ CTT (k) = i[J ffT (0) - J ffT (k)] 



(8) 



To model a simple 1:1 electrolyte one needs only 5 = 2 
components, say + and — with q± = ±q. We may then 
identify basic energy and range scales, jo and i?o, via 



J°_(O) = k B T = 2j >0 t 



(9) 



and, assuming short-range isotropic nonionic couplings, 

A J° + + 2A J% _ + AJ°__ « 2j Rlk 2 > , (10) 

as k —>■ 0. Now, if J++(R) = J (R), the model is 

(fully) charge symmetric (as is the well known contin- 
uum Restricted Primitive Model or RPM [HE A 
suitable charge asymmetry parameter is then 14] 



S , = max lAJa 



.(k)-AJ__(k)|/fc B r , (11) 



which might, e.g., be used to represent distinct ionic sizes, 
a+ and a_ [6(a)]. As the simplest "Basic Ionic Model" it 
suffices to take only nearest neighbor couplings, J+_ > 
and J° + = — J"_; then one has jo = 2 d ~ 1 J°_ and Sj — 
d\J+ + \/2 d ~ 2 \ J+_|. Finally, as a dimensionless measure 
of the relative strength of the Coulomb interactions, it is 
helpful to introduce the ionicity [15j 



lo = q z /a d - 2 k B To 



(12) 



Of course, even this Basic Ionic Model is insoluble for 
d > 2. Accordingly, in the standard way jToL UlL fl2T |. 



we "sphericalize" these multicomponent models by taking 
the spins sf as unbounded continuous variables subject 
only to S spherical constraints, (s 2 ) — 1, enforced with 
the aid of S Lagrange multipliers which, for later conve- 
nience, we write as A CT -I- ij aa (0). In full generality, the 
singular part of the total free energy f[T, h; A(T, h)] in 
the thermodynamic limit, is then [lj| 

/ s (T;A) = ifc B T^ln{|A(k;A)|/(fc B T) 5 } , (13) 

where J k = {a/2-K) d J kgB d d k and h = {h a ) and A = 
(A CT ), while the iSx<S matrix A has elements 

A ffT = [A,, + AJctctW]^ - i(l - S aT )Zr(k) . (14) 

The field-dependent contribution to the free energy is 



fh 



■^hlA-^OjAJIh) with h = 2A(0;A)m, (15) 



and m a 



(.Sct) = 2{n a ) — 1 so that p a 



where m = ( 

|p™ OI (l+m a ). Finally, the Lagrange multipliers A(T, h) 
are determined implicitly via the S spherical constraints 
(s 2 ) = {df /d\ a ) = 1. These results are valid while the 
eigenvalues of A remain positive; the vanishing of any 
one of them signals a thermodynamic singularity. 

For brevity hereon we focus on the two-species, 1:1 
case (with q± = ±q), the + and — ions residing on one 
of two hypercubic sublattices displaced by 4 (a, a, ■ ■ • , a): 
see FigQ]for d = 3. The eigenvalues of A, to be called A z 
and Ajv for reasons soon to be evident, are then simply 
the + and — roots of the quadratic equation |A— acl| = 0, 

A x (k;A) = A + AJ(k)-(-) ,Jx D(k;A) with (16) 
2? 2 (k; A) = [At + Ajt(k)] 2 + i Jf_(k) , (17) 

where, for each variable g+, <?++, etc., we have 
introduced the mean g — \{g+ + g~) and difference 

= j (.9+ — 9-), while the square root is chosen so that 
T>(0;X) > and A x (k) is analytic. 

A pivotal result now transpires ^| , namely, the linear 
decomposition of the structure factors via 



5 xr (k;A) _ B£ y (k;A) , B z XY (k;\) 



k B T/Apa d 



in which Bl 



Bl 



A N (k; A) 
= l-B* 



A z (k; A) 



B N ZZ , while 



K«(k;A) 



4|J+_(k)|/X>(k;A) 



(18) 



(19) 



and B^ z + B^ z = 0. In the charge symmetric cases, 
all the g< variables vanish by definition, so B™ N = 1 
and B^ N — B zz = 0: this implies that S NN is entirely 
governed by A N , and, likewise, S zz by A z , so justifying 
the notation. Conversely, in nemsymmetric models, both 
A N and A z contribute to all the structure functions. 
Finally, the spherical constraints reduce, first, to 

k B TJ d (\) + h 2 /4X 2 = 1 with X — X - j' , (20) 
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where j' Q — (0) ^ij], while the basic integral 

J d (A)-i^[A; 1 (k;A)+A- 1 (k;A)] , (21) 

becomes singular, typically as Aj7d ~ —A 1 / 7 , when A (~ 
V) vanishes, where 7 is the appropriate critical exponent 
0,E1E3; lastly At is given (implicitly) by 

At f 1/ |A(k; A) I = - / Ajt(k)/ |A(k; A)| . (22) 

Recall, however, that At vanishes identically in charge 
symmetric cases: more generally, this result relates A* to 
the asymmetry parameter Sj [see Eq. ifTDl ] via 

Xl-csiJ^S,,, (23) 

where, however, cs[Jar] might vanish "accidentally". 

We are now in a position to answer, with explanations, 
the questions posed after Eqns. (QJ-®- For specific nu- 
merical results we will invoke the Basic Ionic Spheri- 
cal Model (BISM), i.e., the sphericalized version of the 
d = 3, S = 2, 1:1 model with p+ = p_ = \p and nearest- 
neighbor interactions as set out with Eqns. il9")l- H12ll . 
When 8 = Iq = one readily finds from the vanish- 
ing of A N (0) when A — > 0,0, that standard 
spherical-model criticality and scaling pertains, with ex- 
ponents (3 = 5 and 7 = 1 — a = (2 — rf)u, where short- 
range (non-Coulomb) interactions lead to 77 = and 
7 = max{2/(d - 2); 1} (for d > 2). 

For Iq and Sj not too large, the same situation pre- 
vails — contrary to speculations (for hard-core continuum 
ionic models) that lack of symmetry might lead to mean- 
field criticality This follows most directly from the 
small k behavior of the eigenvalues, namely, 

A N (k; A) = A + j R 2 N {\) k 2 + 0{k 4 ) , 

A z (k; A) = (S d q 2 /4a d )[l/k 2 + R% + 0(k 2 )} , (24) 

where S d is the area of a unit d-sphere, while 0] 

R 2 N (X) = R 2 - S' d a 2 l Q - 2a 2 At 2 /S d I j 2 , 
i? 2 (A,k) = 2 a 2 (A + ^)/S d Z j + a 2 £ 4 (k), (25) 

in which S' 3 = 7r/144 and £4 = J2 a ^a/^ 4 - The crucial 
feature, following directly from llKit . is that the Coulomb 
singularity, characterized by ^°(k) ~ 1/fc 2 , cancels out 
of A N exactly thanks to electroneutrality. (The absence 
of this possibility for S — 1 results in the destruction of 
gas-liquid criticality by Coulomb interactions [11(b)].) 

To ensure the stability of the critical point, one also 
needs Q: (i) Iq < If ax , with If ax = ffvr ~ 2.77; 
(ii) R 2 N > 0, which restricts Iq and Sj to the interior 
of an ellipsoid with a vertex at Sj = Iq = 0, which, 
for d = 3, is ±S 2 + [±tt! - (R 2 /a 2 )} 2 < R 4 /a 4 ; 
and (iv) the absence of competing minima in A N (k), 



which is satisfied for sufficiently small Sj, specifically 
by Sj < 1 - j^ttIo in the BISM. The solution At (A) 
of (fSSj) then varies smoothly when A — > and (123)) ap- 
plies. By lj2f)l) criticality is restricted to h = and occurs 
at k B T c = 1/J d {0, Aj) and p c = a~ d = \p max - For the 
BISM, we find T c » T^/Kj, cc , to lowest order in Sj and 
Iq, with Kbcc ~ 1.39 0|- It transpires that T c (Sj) is 
a decreasing function of the asymmetry in accord with re- 
cent simulations of hard-core continuum electrolyte mod- 
els [6(a)] that, however, contradict various approximate 
theories. Fu rthermore, a ter m v arying as Iq 1 (in d = 3) 
appears ^ij], in accord with and in analogy with [l7j . 
As regards the density correlation lengths, i|18ll yields 

^ tl (T,p) = R N [A(T»] [jo/A(T,p)] 1/2 , (26) 

for all acceptable Iq and Sj; when p = p c , this diverges 
as CSr/i". Furthermore, all higher moments [see in- 
cluding the "true" correlation length £ Nj00 satisfy 
Snj>/Z*,i -» 1 when (T,p) -» (T c ,p c ). For the BISM, 
with small Iq and (5j we have t;° N ~ a/irKb cc ~ 0.229 a, 
close to the d = 3, nearest neighbor Ising model value 

m. 

By contrast, the near-critical charge correlation lengths 
depend radically on symmetry. In charge symmetric 
models, where B zz = 1, it follows from II18II and 12411 
that the Stillinger-Lovett sum-rule © is valid for all fluid 
regimes including the critical point. However, the true 
charge screening length is given by 

6.00 (T,p) = R z (T,p)[l + 0(I 2 )} , (27) 

where, from l|25|l we find R z /£,d — » 1 as p — » whereas 
near criticality one has R z ,cHd,c ~ 2\/Tq/T c ; that yields 
R z ,c/in,c ^ 2.36 for the BISM (for Iq and Sj not too 
large). Furthermore, when t — ► 0, the screening length, 
£z,oo{pc,T) gains, in general, a singular correction factor 
[1 + ci-ai 1 -"] Up to 0{Tq), the higher moments of 
S zz satisfy (£ Z , P ) P ~ ioR^" 1 ■ 

On the other hand, in nonsymmetric cases B zz in l|18|) 
does not vanish : rather one has 

B~ z (k; A) = 4At 2 fcV/<? 2 Z 2 j 2 [l + 0(k 2 )] > . (28) 

Consequently, all charge correlations become infected by 
the divergent density fluctuations controlled by A N (k). 
Nevertheless, because of the factor k 4 , the Stillinger- 
Lovett relation © remains valid in the fluid re gim e ex- 
cept at criticality where it fails and we find 0Uil 

(k,i/&)c = 1 + - 1 + w 2 c cjS 2 + ..., (29) 

where w 2 = 2a 2 / Sd1oR 2 N j 2 , and, recalling <2.H . we note 
that C5[J CTT ] 7^ for the BISM. This critical point failure 
implies a breakdown of full screening that is necessarily 
associated with slow decay of certain ionic correlations 
0. Indeed, when (T,p) — > (T c ,p c ) the charge decay 
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TABLE I: Charge correlation lengths near criticality where 
£jv diverges, while R z /£,d = 0(1). The ionicity, T oc q 2 , and 
asymmetry factor w c csSj, are defined in 11211 . 12311 and 12911 . 



6m 



SZ,p 



charge symmetric 
: £ D = (47V P q 2 /k B T) 1 / 2 , 
the Debye length, 

e D R%{T,p)=o{&), 
aRi^ii+oiio)], 



non-symmetric 
= £n for {T,p)^(T c ,p c ), 
> Cd,c at (To, pc), 
= -&l(v) c cs8j)*& - 



-(w c cs5„ 



\2 e 2 t2(p-l) 



fl z [l+O(T 2 )] = Ofo), « C« ~ 1/f 



length £z !00 asymptotically approaches the density cor- 
relation length, £n,oo, and thus diverges as (for 
P = Pc)- 

However, the fourth charge correlation moment is given 
(except at criticality) by 



(30) 



so that |£ z . 2 (/O c , T)\ diverges more weakly as 1/t"/ 2 . Note 
also that on approaching criticality, £ z 2 changes sign 
(with respect to the symmetric case); for small 8,, this 
crossover occurs in the BISM at t = t x « 0.265 <5, 7 . 
More generally the higher moments in nonsymmetric sys- 
tems satisfy £ z p p w — $$f<Z c w c^c leading to a hier- 
archy of critical exponents ~ \/t v V-~ 1 / p ' . Notwith- 
standing the divergence of the charge correlation length, 
iz,oo ~ £w ~~ > °° ; the charge-charge pair correlation func- 
tion, G zz (r), decays exponentially at (T c ,p c )\ Indeed on 
approach to criticality we obtain 



G zz (r) cc2# 



2 To c§ e" 



T c J 2 i? 2 &(2»7 



Rt r 



(31) 



to leading orders for d = 3. At T c only the second term 
survives since the first vanishes as t Av . 

Finally, as regards the cross-correlations embodied 
in the charge-density structure factor S NZ , we find 
S NZ (k — 0; T, p) = except at (T c ,p c ) where the value 

2 

(At ^o/jo-R«)c is realized. Moreover, on defining mo- 
ments in analogy to (JHJ , one obtains 

C P Z , P (T, p) « (\^l/3oRl)e N P oc(T, P) , (32) 

near criticality; of course, all these moments vanish iden- 
tically in charge symmetric models since = 0. 

In summary, we have analyzed a class of exactly solu- 
ble spherical models for 1:1 ionic systems and shown that 
the Coulomb interactions do not change the gas-liquid 
critical universality class — contrary to some suggestions 
[HQ- The couplings between the charge correlations and 
the divergent critical density fluctuations 0, 0] follow 
mainly from a remarkable structure-function decomposi- 
tion, Eq. ltl8|l . that respects the Stillinger-Lovett (SL) 



sum rule (unlike 5]). Our principal results are collected 
in Table U : they are broadly consistent with Ornstein- 
Zernike-based arguments advanced for hard-core contin- 
uum electrolytes [3] . In charge-symmetric models density 
fluctuations are not directly coupled to two-point charge 
correlations which remain of short-range and obey SL 
near and at criticality; but in more realistic nonsymmet- 
ric systems the density fluctuations "infect" the charge 
correlations which hence exhibit the same diverging cor- 
relation length. Moreover, the SL rule is then violated 
at criticality fl4lllfi| indicating an anomalous conducting 
state IsB. 
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